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Abstract 



A study of charm fragmentation into D*^ and Df in e^e annihilations at 
= 10.5 GeV is presented. This study using 4.72±0.05 fb^^ of CLEO II data 
reports measurements of the cross-sections a{D*^) and a{Df) in momentum 
regions above x = 0.44, where x is the Dg momentum divided by the maximum 
kinematically allowed Dg momentum. The Dg vector to vector plus pseu- 
doscalar production ratio is measured to be Pvix{Df) > 0.44) = 0.44 it 0.04. 
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I. INTRODUCTION 



The production cross-sections of qq pairs in e^e~ annihilations can be calculated using 
QCD, but the process of fragmentation whereby hadrons are formed is non-perturbative and 
phenomenological models are used to describe it. Two properties of hadron production that 
can be experimentally measured are the hadron momentum distribution and the relative 
population of available spin states. 

Measurements of primary hadron fragmentation can be challenging due to cascades from 
higher order resonances that can be indistinguishable from the primary hadrons. The study 
of and Z)*"*" fragmentation in e^e~ annihilations at ^/s = 10.5 GeV benefits from the 
fact that L = 1 charm mesons have not been observed to decay to either or D*^ 
and the influence of B events is kinematically eliminated for x{Dg) > 0.4, where x is the Ds 
momentum divided by the maximum kinematically allowed Dg momentum. The Dg system 
is thus particularly well suited for the measurement of the vector to pseudoscalar production 
ratio. 

The vector to pseudoscalar production ratio is usually described using the variable 

where P and V represent, respectively, the number of pseudoscalar and vector mesons directly 
produced through a particular production mechanism, e.g. e~^e~ annihilations. Counting the 
number of spin states available to an L = meson leads to the expectation that Py = 0.75. 
This spin counting model has been shown to be useful for describing the D*^ spin align- 
ment [0, but most measured values of Py have been significantly lower than 0.75 for charm 
mesons. Other models based upon the mass difference between the vector and pseudoscalar 
states predict values of Py that are less than 0.75 0, but more precise measurements are 
needed to better determine any relationship between Py and the mass difference. 



II. DETECTOR AND EVENT SELECTION 

The data in this analysis were collected from e+e^ collisions at the Cornell Electron 
Storage Ring (CESR) by the CLEO II detector. The CLEO II detector is a general purpose 
charged and neutral particle spectrometer described in detail elsewhere The dataset 
used in this analysis contains 3.11 ± 0.03 fb~^ of data collected at the T(4S') resonance and 
1.61 ± 0.02 fb-i of data collected below the bb threshold (about 60 MeV below the T(4S') 
resonance), for an approximate total of 5 x 10^ cc events. 

In this analysis, D*'^ mesons are reconstructed via the decay D*^ —>■ D^'j and D'^ 
mesons are reconstructed via the decay chain D'^ — > 07r+ with — >■ K^K^ (inclusion of 
charge conjugate modes is implied throughout this paper). 

All charged tracks used in this analysis are required to have an origin close to the e^e~ in- 
teraction region and must be well reconstructed. When drift chamber particle identification 
information is available, the specific ionization, dE/dx, must be within two standard devia- 
tions of the expected value for candidate kaon tracks and within three standard deviations 
of the expected value for candidate pion tracks. 
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Showers in the crystal calorimeter are considered as photon candidates if they have a 
minimum energy of 100 MeV, are within either the barrel (| cos^s| < 0.71, where Og is the 
angle between the shower and the e"*" beam direction) or endcap (0.85 < |cos^<j| < 0.95) 
regions, have an energy deposition consistent with that expected for a photon, and do not 
include any crystals near a projected charged track. 

Candidate mesons are reconstructed using all appropriately signed combinations of 
candidate kaon tracks in an event. The invariant mass M{KK) is required to be within 
8.4 MeV/c^ (approximately 2 standard deviations) of the known mass |l[]. Candidate 
mesons are reconstructed using all combinations of candidates and candidate pion 
tracks in an event. Candidate D*~^ mesons are reconstructed using candidate photons, and 
(f)7T combinations with invariant mass M{KK7c) within 20 MeV/c^ (approximately 2.5 to 3 
standard deviations) of the known mass. 

Because the must be polarized in the helicity-zero state ina Dg 0vr decay, the decay of 
the has an angular distribution proportional to cos'^a, where a is the angle between the 
and -D^ momentum vectors in the rest frame. Since the background angular distribution 
is flat, the signal to background ratio is improved by requiring |cosq;| > 0.35. The signal to 
background ratio is further enhanced by requiring that cos^^ > —0.8, where 6*^ is the angle 
of the TT momentum vector in the rest frame relative to the momentum vector in the 
laboratory frame; the signal distribution is flat in this variable while background events peak 
at cos 6'^ = —1.0. Because of the minimum energy restriction for photon candidates, signal 
photons traveling in a direction opposite to the D*'^ direction in the laboratory frame are 
excluded from the candidate sample. By requiring cos6'^ > —0.8, where 9^ is defined as the 
angle of the photon momentum vector in the D*'^ rest frame relative to the Z)*^ momentum 
vector in the laboratory frame, additional background Z)*"*" candidates are suppressed. 

Low momentum candidates are difficult to analyze because of the large amount of 
background from combinatorics as well as B decays. The analysis is therefore restricted to 
x{Df) > 0.44 where 

and 



Pmax{,Dt) = ^El^^-ml+ . (3) 
For D*~^ candidates, the x{D^) requirement is replaced by x{Dl^) > 0.5 where 

fmax \-^s I 

and 

Vmax 



\ 



1 1 

ml* I — m 



2 



Ebeam-^ ^1 -ml,. (5) 



In principle, B -D^vr can result in xiD^) ~ 0.5. However, such decays are b ^ u 
transitions and thus heavily suppressed, so they are expected to be a negligible source of 
background. 
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Based on the assumption that all observed D*~^ are primary, the D*~^ momentum spec- 
trum is simply studied by measuring the D*~^ yield in eight equal sized bins of x{D*~^) over 
the range 0.5 < x{D*~^) < 0.98. However, the observed can be primary or D*~^ daugh- 
ters. In order to study the momentum distribution of primary mesons, it is necessary 
to subtract out the D*~^ contribution to the yields. Since all D*~^ are assumed to decay 
to D'^, the yields from D*~^ decays can be accounted for by simply measuring the D*'^ 
yields as above, but in bins of the variable x{Df) rather than x{D*~^). After the D*'^ yields 
are corrected for efficiency and the branching ratio B{D*~^ they are subtracted 

from the efficiency corrected yield in each x{D^) bin to calculate the primary yield. 

III. FITTING 

The yields are projected onto AM = M^KKti-i) - M{KKtt) for D*-^ candidates 
and the D'^ yields are projected onto M{KK7c) for candidates. Fitting shapes for the 
peaks in these distributions are determined using a sample of Monte Carlo events generated 
using the Lund jetset 7.3 program combined with a GEANT-based CLEO II detector 
simulation, where every event contains a D*"*" or decaying through the modes specified 
above. 

The AM distributions in data and the signal Monte Carlo sample are simultaneously fit 
to the sum of an asymmetric Gaussian for the signal and separate second-order Chebyshev 
polynomials for the background in each distribution. An asymmetric Gaussian is used be- 
cause of the larger tail on the lower side of the peak attributable to energy leakage in the 
calorimeter. The fits to data used to determine the D*^ yields in the selected regions of 
x{D*~^) and x{Df) are shown in Figs. |l] and 

The M{KKti) distributions in data and the signal Monte Carlo sample are simultaneously 
fit to the sum of a double Gaussian with common mean for the signal, a Gaussian for the 
signal, two straight lines joined by a quadratic for the combinatoric background in data, 
and a first order Chebyshev polynomial for the small amount of background in the Monte 
Carlo sample. The fits to data used to determine the yields in the selected regions of 
x{D'^) are shown in Fig. ^ 

IV. EFFICIENCIES 

The and D*+ detection efficiencies are estimated using a sample of Monte Carlo 
events that contains signal as well as background events and is independent of the signal 
Monte Carlo sample used in the fitting procedure. The D^^ efficiency values in the x{D*'^) 
regions are listed in Table |, while the and D*"*" efficiencies in the x{D~^) regions are 
listed in Table 0. 

For the D*"'" production study, the efficiency for each x(Z^*^) bin is measured using 
the fitting procedure described above. The binned raw efficiency values within the range 
0.50 < x{Dl^) < 0.98 are fit with a first order Chebyshev polynomial to provide a smoothly 
varying efficiency as a function of x{D*~^). The smoothed efficiency value at the center of 
each x{D*~^) region is used to calculate the efficiency corrected D*~^ yield and cross-section. 
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FIG. 1. Fits to the AM = M{KK-K-j) - M (KKtt) distributions for candidate D*+ events that 
are used to determine the D*^ yields in the eight x{D*'^) ranges (a) 0.50 — 0.56, (b) 0.56 — 0.62, 
(c) 0.62 - 0.68, (d) 0.68 - 0.74, (e) 0.74 - 0.80, (f) 0.80 - 0.86, (g) 0.86 - 0.92, and (h) 0.92 - 0.98. 



7 



80 



3110300-002 





40 









80 




An 
40 


o 

















in 






40 


(A 




c 




> 





LU 






20 









4 





1 


1 1 




— 1 

\a 


1 \— 

iJ^ * - 




'''' +*f^^**t^ 


-"/ 


\ 


1 


(g) L 


[ : 


— ^ — r 
(i) 


1 

m 1 


— 


T TT Tm T T T TTTTT 

1 1 1 


1 



60 
40 

20 


60 

40 

20 



40 

20 



10 



4 






(j) 



0.10 0.14 0.18 0.22 0.10 0.14 0.18 0.22 
AM = M(KK7r7)-M(KK7r) (GeV/c^) 

FIG. 2. Fits to the AM = M{KKttj) - M{KK-k) distributions for candidate D*+ events 
that are used to determine the D*+ yields in the ten x{Df) ranges (a) 0.44 — 0.50, (b) 0.50 — 0.56, 
(c) 0.56 - 0.62, (d) 0.62 - 0.68, (e) 0.68 - 0.74, (f) 0.74 - 0.80, (g) 0.80 - 0.86, (h) 0.86 - 0.92, (i) 
0.92 - 0.98, and (j) 0.92 - 1.00. 
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FIG. 3. Fits to the M{KK7r) distributions for candidate Df events that are used to determine 
the D+ yields in the ten x{Df) ranges (a) 0.44-0.50, (b) 0.50-0.56, (c) 0.56-0.62, (d) 0.62-0.68, 
(e) 0.68 - 0.74, (f) 0.74 - 0.80, (g) 0.80 - 0.86, (h) 0.86 - 0.92, (i) 0.92 - 0.98, and (j) 0.92 - 1.00. 
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x{D*~^) region D*+ Efficiency 

0.50 - 0.56 0.172 ±0.008 

0.56 - 0.62 0.175 ± 0.006 

0.62 - 0.68 0.178 ± 0.005 

0.68 - 0.74 0.181 ± 0.004 

0.74 - 0.80 0.184 ± 0.004 

0.80 - 0.86 0.187 ±0.005 

0.86 - 0.92 0.190 ±0.006 

0.92 - 0.98 0.193 ±0.008 

(0.92 - 1.00) (0.194 ± 0.008) 

TABLE I. detection efficiencies in the specified regions of x{D*'^), where the efficiency has 
been smoothed using a fit to the raw efficiency spectrum with a first order Chebyshev polynomial. 

The values that are in parentheses are used only for the calculation of the total yield and 

cross-section for x(D*+) > 0.5. 



x{Df) region 


Df Efficiency 




D*+ Efficiency 


0.44 - 0.50 


0.366 ± 0.010 




0.180 ± 0.009 


0.50 - 0.56 


0.369 ± 0.008 




0.179 ± 0.007 


0.56 - 0.62 


0.373 ± 0.007 




0.178 ± 0.006 


0.62 - 0.68 


0.376 ± 0.005 




0.177 ± 0.005 


0.68 - 0.74 


0.379 ± 0.005 




0.176 ± 0.004 


0.74 - 0.80 


0.382 ± 0.005 




0.175 ± 0.005 


0.80 - 0.86 


0.386 ± 0.007 




0.174 ±0.006 


0.86 - 0.92 


0.389 ± 0.008 




0.140 ±0.027 


0.92 - 0.98 


0.392 ± 0.010 




0.248 ±0.107 


(0.92 - 1.00) 


(0.393 ± 0.011) 




(0.248 ± 0.106) 


TABLE II. Df and D*^ detection efficiencies in the specified regions of x{Ds) where the 
efficiency has been smoothed using a fit to the raw efficiency spectrum with a first order Chebyshev 
polynomial. The values that are in parentheses are used for the calculation of Pv{Ds). The -D*"*" 
efficiency values with x{Df) > 0.86 are excluded from the smoothing process because they are not 
expected to be modeled by the same function used for x{Df) < 0.86. 
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Measured 


Efficiency Corrected 




x[D* ) region 


JJg 1 leiu 


I leid. 




U.oU - U.OD 


OA 1 1 O 

84 ± 18 


/loo 1 1 1 

488 ±116 


O.IU ± U.Uz 


U.OD - U.o2 


159 ± 19 


ni 1 1 1 o o 

911 ± lo8 


0.19 ± U.U3 


O.bz - U.oo 


176 ± 19 


972 ± 137 


0.21 ± 0.03 


n Ra n 7/i 

U.Do - U. /4 


1 4-17 


1 ni /I 4- 1 Qi 


n 91 4- n IT? 


0.74 - 0.80 


202 ± 17 


1095 ± 135 


0.23 ± 0.03 


0.80 - 0.86 


172 ± 16 


918 ±117 


0.19 ±0.02 


0.86 - 0.92 


148 ± 13 


778 ± 100 


0.16 ±0.02 


0.92 - 0.98 


84± 10 


434 ± 65 


0.09 ±0.01 


(0.92 - 1.00) 


(99± 11) 


(511 ± 74) 


(0.11 ±0.02) 


0.50 - 1.00 


1219 ± 47 


6687 ± 260 ± 497 


1.42 ±0.06 ±0.11 


TABLE III. D*^ 


yields, efficiency corrected yields and cross sections 


in the specified regions 


of x{D*~^), where B 




^ (j)TT+)B{(j) K+K-). 


The uncertainty in the 



efficiency corrected yields and cross-sections includes both statistical and systematic error. When 
two errors are presented, the first is statistical while the second is systematic. 



For the fragmentation study, the and D*^ efficiencies are measured in each 
x{D^) bin using the fitting procedure described above. The binned raw efficiencies 
within the range 0.44 < x{Df) < 0.98 are fit with a first-order Chebyshev and the smoothed 
efficiency values are used to calculate the efficiency corrected yields. The binned raw 
D*~^ efficiencies with 0.44 < x{Df) < 0.86 are fit with a first-order Chebyshev polynomial 
but the efficiencies in the region x{Df) > 0.86 are excluded from the fit because of expected 
efficiency loss due to the larger proportion of photons in that region with energies less than 
100 MeV. The smoothed efficiency values are used to calculate the efficiency corrected D*~^ 
yields for x{Df) < 0.86, while the raw efficiency values are used for x{Df) > 0.86. 



V. RESULTS 

The D*~^ yields, efficiency corrected yields and cross-sections in the eight x{D*^) regions 
are all listed in Table |T|. These same quantities for and D*'^ in the nine x{D^) regions 
are listed in Tables |^ and 0, respectively. The calculated primary yields and cross- 
sections are presented in Table 0. 

By summing the efficiency corrected D*+ and primary yields listed in Tables fTT| and 
f7\ respectively, Eq. ([l|) could be used to calculate Py for x{D'f'>^) > 0.5. However, the 
uncertainties in the D*'^ yields are essentially counted twice due to the subtraction used to 
calculate the primary yields. Py can however be calculated in a way that avoids this 
subtraction. Since all observed D*^ mesons are assumed to be primary, all observed Z)+ 
mesons are assumed to be either primary or D*'^ daughters, and all D*^ are expected to 
decay to a D'^, Eq. (P can be rewritten as 

P -11EI11 IR\ 

Pv-Y{Di)' "'^ 
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Measured 


Efficiency Corrected 




x{Df) region 


D+ Yield 


D+ Yield 


B-a{Dt){vh) 


0.44 - 0.50 


546 ± 43 


1491 ± 140 


0.32 ± 0.03 


0.50 - 0.56 


663 ± 40 


1795 ± 139 


0.38 ± 0.03 


0.56 - 0.62 


933 ± 40 


2504 ± 160 


0.53 ±0.03 


0.62 - 0.68 


1019 ± 38 


2712 ± 160 


0.57 ±0.03 


0.68 - 0.74 


1080 lb 40 


2847 lb 166 


0.60 ± 0.04 


U. / 4 - U.oU 




0/11 Q 4_ 1 /I K 


U.Oi It U.Uo 


U.oU - U.OU 




1968 ± 122 


0.42 ± 0.03 


0.86 - 0.92 


404 ± 24 


1038 ± 79 


0.22 ± 0.02 


0.92 - 0.98 


170 ± 14 


433 ± 42 


0.09 ± 0.01 


(0.92 - 1.00) 


(187 ± 17) 


(476 ± 50) 


(0.10 ±0.01) 


0.44 - 1.00 


6516 ± 106 


17250 lb 281 lb 528 


3.65 ±0.06 ±0.11 


TABLE IV. 


Df yields and 


cross-sections in the specified regions 


of x{Df), where 


B = B{Dt ^ <? 


^7r+)B{(t) K+K-). 


, The yields are efficiency corrected using smoothed ef- 


ficiency values. The uncertainty in the efficiency corrected yields and cross-sections includes both 



statistical and systematic error. When two errors are presented, the first is statistical while the 



second is systematic. 




Measured 


Efficiency Corrected 




,x(D+) region 


D*+ Yield 


D*- Yield 




0.44 - 0.50 


104 ± 21 


577 ± 129 


0.13 ±0.03 


0.50 - 0.56 


145 ± 19 


808 ± 133 


0.18 ±0.03 


0.56 - 0.62 


200 ± 20 


1122 ± 152 


0.25 ± 0.03 


0.62 - 0.68 


181 ± 18 


1020 ± 138 


0.23 ± 0.03 


0.68 - 0.74 


208 ± 18 


1180 ± 148 


0.27 ± 0.03 


0.74 - 0.80 


182 ± 16 


1040 ± 133 


0.23 ± 0.03 


0.80 - 0.86 


149 ± 14 


852 ± 113 


0.19 ±0.03 


0.86 - 0.92 


70 ±9 


501 ± 202 


0.11 ±0.05 


0.92 - 0.98 


16 ±5 


63 ±39 


0.01 ±0.01 


(0.92 - 1.00) 


(15 ±5) 


(60 ± 37) 


(0.01 ±0.01) 


0.11 - 1.00 


1253 ± 49 


7160 ± 279 ± 550 


1.61 ±0.06 ±0.11 



TABLE V. D*'^ yields and cross-sections in the specified regions of x{Df), where 
B = B{Df (f)Tr'^)B{^ K^K~). The first seven D*+ yields are efficiency corrected using 
the smoothed efficiency values while the efficiencies for x{Df ) > 0.86 arc corrected using the raw 
efficiency values. The uncertainty in the efficiency corrected yields and cross-sections includes both 
statistical and systematic error. When two errors are presented the first is statistical while the 
second systematic. 
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Jrilllicxiy J-^ g 1 IciQ 
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0.74 - 0.80 


1315 lb 204 


0.28 ± 0.04 


0.80 - 0.86 


1063 ± 173 


0.23 ±0.04 


0.86 - 0.92 


506 ± 229 


0.11 ±0.05 


0.92 - 0.98 


370 ± 57 


0.08 ±0.01 


(0.92 - 1.00) 


(417 ± 64) 


(0.09 ±0.01) 


0.44 - 1.00 


9652 ± 408 ± 760 


2.05 ±0.09 ±0.16 



TABLE VI. Calculated primary D'^ yields and cross sections in the specified regions of x{D^). 
The cross-section is presented as iS-u where B = B{Df — > (f)Ti^)B{(j) K^K^). The errors in each 
x{D'^) bin include both statistical and systematic uncertainty. When two errors are presented, the 
first is statistical while the second is systematic. 



where T(M) is the total number of M mesons in the CLEO II data sample. In terms of the 
quantities measured using the decay modes chosen for this analysis, 



where n{M) is the efficiency corrected yield of M mesons in a particular x{D^) region. 
Using this method, Pv{,x{D+) > OAA)B{D*+ Df-f) = 0.42 ± 0.02. Using the value 
B{D*+ Dti) = (94.2±2.5)% § leads to Pv'(a;(£'s) > 0.44) = 0.44±0.02(stal)±0.01(6r.). 



VI. SYSTEMATIC UNCERTAINTY 

The systematic error for the total D'^ and D*"^ yields is determined by varying the 
selection and fitting procedures as described below and taking the variance in the total yield 
as the estimate of the error. The variance is also determined on a bin-by-bin basis and 
the average percentage variance in the individual bins is taken as the estimated systematic 
uncertainty for all bins. The uncertainties in the D*'^ yields for the range 0.86 < x{D'^) < 1.0 
are averaged separately since those values are not smoothed and the errors are quite large 
due to the limited number of D*^ events in that region. Systematic uncertainties on the 
various yields are listed in Tables |V1I| and |VIl.i| . 

The acceptance angles for showers implicitly alter the acceptance of tracks since there 
is a high degree of correlation between the fiight directions of the and the photon in 
the detector. There is also a correlation between the photon energy and the decay angle 
of the D*'^. Varying the shower acceptance angles to | cos 6*3! < 0.5 changes the total D*^ 
yields and the bin-by-bin yields by approximately 6%, while changing the minimum shower 
energy to either 90 MeV or 110 MeV changes the total D*~^ yield by approximately 2% 
and the bin-by-bin yields by approximately 3%. A 3% overall systematic uncertainty in 
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Variation 


Percent Variance 


peak in M{KKtt) fit with Gaussian 


2% (3%) 


peak in M^KKtt) fit with double Gaussian 


1 0/ /1 0/ \ 


background fit with quadratic 


2% (3%) 


TABLE VII. Percent variance in total(bin-by-bin) yields 


compared to the nominal yields 


due to the hsted sources of systematic error. 




Variation 


Percent Variance 


15 MeV/c^ wide M{KKtt) signal region 


1%(3%) 


25 MeV/c^ wide M{KKtt) signal region 


1%(2%) 


AM peak fit with Gaussian 


1%(1%) 


AM peak fit with double bifurcated Gaussian 


2%(3%) 


cos 6s < 0.6 


6% (6%) 


E{j) > 90, 110 MeV 


2% (3%) 


Uncertainty in 7 efficiency from 13{r] ^^)/B{ri — > Stt*^) study 


3%(3%) 


TABLE VIII. Percent variance in total(bin-by-bin) D*^^ yields 


compared to the nominal yields 


due to the listed sources of systematic error. 





photon reconstruction has been estimated by comparing the world average value of B{ri —>■ 
77)/i3(?7 — 37r°) |l|] with the relative yields of — > 77 and rj Sn^ in data and Monte Carlo. 

Additional uncertainty exists because of differences in invariant mass distributions be- 
tween data and Monte Carlo and possible inadequacies of the fitting functions used to deter- 
mine the yields. This uncertainty is estimated by altering the fitting shapes used to obtain 
the and D*~^ yields. Varying the fitting technique for the M{KK7c) projections by e.g. 
using a Gaussian for the signal peak, a double Gaussian with common mean for the 
signal peak, or a second-order polynomial for the background alters the total Df yield by 
approximately 3% and the bin- by-bin yields by approximately 4%. Using a single Gaussian 
or double bifurcated Gaussian with a common mean for the peak in the AM distribution 
alters the total D*~^ yield by approximately 2% and the bin-by-bin yields by approximately 
3%. 

There is also an uncertainty related to the requirement that M{KK7r) be within 20 
MeV/c^ of its nominal value. Widening this requirement to 25 MeV/c^ and narrowing 
it to 15 MeV/c^ has resulted in an approximate 2% error in the total D*^ yield and an 
approximate 4% error in the bin-by-bin yields. 

The uncertainties in the efficiency values shown in Tables | and |T| vary for each region 
of x{Df) and x{D*~^) due to limited Monte Carlo statistics and the smoothing process. 
For instance, the errors in the smoothed efficiency values near the limits of the x region 
studied are higher than those in the middle of the region due to the uncertainty in the slope 
of the function used in the smoothing process. The errors in the efficiency contribute to 
the systematic uncertainty on a bin-by-bin basis and the percentage errors are added in 
quadrature for the determination of the percentage error for the total yields. 

All of the individual systematic uncertainties associated with a given yield are added 
together in quadrature with the percentage error in the efficiency to determine the total 
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Fit Results X^/d.o.f. 



Andersson: 

D*+: a = 0.9±0.2, m_L = 1.7±0.1 1.9/5 

D+: a = 1.1±0.2, m_L = l.SibO.l 3.2/6 
Peterson: 

D*+: ep = 0.056 ±0.008 20.5/6 

D+: ep = 0.10 ±0.02 17.4/7 

TABLE IX. Results of fits to the D*'^ and Df spectra in x{D*^) and x{D^), respectively, 
with the Anderson et al. and Peterson et al. analytical fragmentation functions. 

systematic uncertainty in the D'^ and D*'^ yields. These systematic uncertainties are already 
included in the errors in the yields in Tables |ITT|, ^V\, |V] and |V1[ After including the total 
systematic uncertainty, Pvix{Df) > 0.44) = 0.44 ± 0.02{stat.) ± Om{syst.) ± 0.01(6r.). 



VII. DISCUSSION OF RESULTS 

The momentum distributions of hadrons created in the fragmentation process are com- 
monly modeled with either the Andersson et al. symmetric fragmentation function [0 or the 
Peterson et al. fragmentation function [If. Both of these functions depend upon z = , 
where Eh is the energy of the hadron, is the hadron momentum parallel to p, the momen- 
tum of the primary quark from the production process, and E is the energy of the primary 
quark. The Andersson function is 

f{z) oc z-\l - zf exp{-b ml/z) , (8) 



where a and b are free parameters, m± = \Jm'^ + P±i is the mass of the primary quark 
and p_[_ is the hadron momentum perpendicular to p. The Peterson function is 

z[l - [l/z) - ep/[l- z)f 

where ep is the single free parameter. 

To properly compare fragmentation models with data it is necessary to use the above 
functions in a full Monte Carlo simulation that incorporates photon radiation, gluon radiation 
and other effects. To facilitate comparison with other experimental results, x is used as an 
approximation of z and a binned fit to the data is performed using these two functions as 
shown in Figs. ^ and |^. Since the parameters b and m_L only appear in Eq. as a product, 
the constraint 6=1 has been used for the fit, thereby changing the interpretation of the 



value of m^. The numerical results from the fits are listed in Table |^ The normalizations 
of these fits are not used to calculate a value of Py due to differences between x and z that 
are non-negligible in the low Dg momentum regime. 

The fragmentation spectra for charm mesons has been studied previously by the CLEO 
collaboration |^ and input parameters for the Andersson et al. model were determined 
using measured fragmentation distributions for D*"*", D"*", Dg and Ac. A comparison of 
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FIG. 4. B • a{D*^) spectrum fit with the Andersson et al. and Peterson et al. fragmentation 
functions, where B = B{D*+ Dt-f)B{Dt (Ptt+)B{^ K+R-). 
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the data presented here with a Monte Carlo distribution using the parameters determined 
in that study, a = 0.60 and b = 0.52, is shown in Figure The use of these values as 
input parameters for charm fragmentation into Dg mesons clearly result in a momentum 
distribution that is too soft, which is not surprising since P-wave charm meson decays to 
D*^, and were not excluded in the prior study. 

High levels of combinatoric background at low values of x{D^) prohibit a good measure- 
ment of Py for the full range of allowed momenta. Based on Monte Carlo simulations 
and the data presented, approximately 75 — 85% of all D*~^ and D'^ are expected to have 
x{Df) > 0.44 and the value of Py presented here is not expected to differ much from Py(all 
x). 

It is possible to make a model-dependent extrapolation of Pv'(all x) using 



> 0) = . / , (10) 



Qv 



where Qy is the percentage of D*"*" that decay to a with x{Df) > 0.44 and Qp is 
the fraction of primary that have x{D~^) > 0.44. Since only about one fifth of either 
fragmentation spectra lies below x{D~^) = 0.44, and because both distributions approach 
zero smoothly as x{Df) —> 0, the ratio Qy/Qp is expected to be close to unity and to only 
depend weakly upon the chosen fragmentation parameters. 

Using the Andersson et al. model with the parameters a = 0.60 and b = 0.52 results in 
Qv = 0.773, Qp = 0.795, Qv/Qp = 0.972, and from Eq. Py{a\\ x{Df)) = 0.45 ± 0.05. 
Changing the input parameters to provide a harder spectrum has a very small effect on 
Qv/Qp- A distribution created with a = 0.4 and b = 0.9, for example, provides a much 
improved representation of the data and results in Qy = 0.860, Qp = 0.875, Qy/Qp = 
0.983 and Pv'(all x{D^)) = 0.45 ± 0.05. This clearly shows that the dependence of the Py 
extrapolation on the choice of fragmentation parameters is indeed weak. 

Based on the results of varying the input parameters for the two models, a systematic 
uncertainty of 3% is estimated for the model-dependent extrapolation resulting in a final 
extrapolated value of Py{all x{Df)) = 0.45 ± 0.05, which is significantly different than the 
expected result based on spin counting. 

Other measurements of Py for charm and bottom mesons have been presented |TD|-|T?|, 
but it is difficult to make direct comparisons between those results and the one presented 
here because of differences in methodology and center-of-mass energies in the other analyses. 
Nonetheless, measurements of Py{B) are generally close to the spin-counting expectation 
while measurements of Py{D) are well below that value as shown in Table |^. 



VIII. CONCLUSION 

In summary, studies of D*'^ and fragmentation in e'^e" annihilations at ^/s = 10.5 
GeV have been presented. Py{x{Ds) > 0.44) has been measured to be 0.44 ± 0.02{stat.) ± 
0.03{syst.) ± 0.01(6r.). When extrapolated to the entire available momentum region this 
measurement deviates significantly from Py = 0.75, the expected result based on simple 
spin counting. 
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FIG. 6. B ■ (j{D*^) and B ■ (T[Df) for primary Dg compared to the Monte Carlo production 
spectrum using the Andersson et al. fragmentation function parameters a = 0.60 and 6 = 0.52 [9], 
where B = B{D+ (^7r+)e(0 ^ K+R-). 



CoUaboration 



Result 



Pv{Ds) = 0.60 ±0.19 
PviD) = 0.60 ±0.05 
PviD) = 0.57 ±0.06 
PviD) = 0.57 ±0.07 
Pv{B) = 0.76 ±0.10 
Pv{B) = 0.76 ±0.09 



ALEPH PI 



ALEPH lOl 



OPAL |11| 
SLD 
L3 
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TABLE X. Results of previous measurements of Py for heavy quark mesons at other experi- 
ments. All of these measurements used data samples collected at resonance. 
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